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ABSTRACT 

A  mathematical  model  was  developed  to  include  line 
effects  in  the  submarine  hydraulic  system  dynamic 
performance  analysis.   The  project  was  undertaken  in  an 
effort  to  demonstrate  the  necessity  of  coupling  the  entire 
hydraulic  power  network  for  an  accurate  analysis  of  any  of 
the  subsystems  rather  than  the  current  practice  of  treating 
a  component  loop  as  an  isolated  system.   It  was  intended 
that  the  line  model  could  be  coupled  with  individual  elements 
in  modular  computer  programs  to  facilitate  design  and  alter- 
ation decisions. 

The  line  model  was  developed  and  coupled  with  a  submarine 
driving-system  servo- loop  in  order  to  demonstrate  the  above 
premises.   However,  due  to  programming  difficulties 
additional  work  is  necessary  to  successfully  complete  the 
analysis. 
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ABBREVIATIONS 


A  area  (in  ) 


BV  blocking  valve 

c  or  c,  constant 

Cj  orific  flow  discharge  coefficient 

C  ,  coefficient  of  external  leakage  (in  /sec/psi) 

C.,  coefficient  of  cross-port  leakage  (in  /sec/psi) 

cs  centstokes 

d  differential  operator 

D  diameter 

d-,  land  diameter  of  valve  (in) 

D  piston  tail  rod  diameter  (in) 

DV  deceleration  valve 

e  exponential  (2.71828...) 

Fj.  load  force 

o„  degrees  Fahrenheit  -  temperature  (lb) 

K  loss  coefficient  due  to  abrupt  area  change 

K  null  flow  pressure  coefficient  (inVsec/psi) 

L  length  (in) 

L  total  cylinder  length  (in) 

Le  effective  length  (in) 

L  length  available  for  ran  to  travel  with  ram 


o 


in  center  position  (in) 


m  .,  mass  of  oil  (slugs) 

m  mass  of  ram  (slugs) 

MT  total  mass  (slugs) 

MIL-L  military  lubricant  specifications 

P  pressure  (psi) 

• 

P  =  dP/dt  change  in  pressure  with  respect  to  time  (psi/sec) 

P  reference  pressure  (psi) 

psi  pounds  per  square  inch  -  pressure 

AP  differential  pressure  (psi) 


~3 

Q  ,  external  leakage  flow  (In  /sec) 

Q--,  internal  cross  port  flow  (in  /sec) 

Q.  supply  flow  from  header  (in~Ysec) 
Q 


return  flow  to  header  (in  /sec) 
SAE  Society  of  Automotive  Engineers 


SG  specific  gravity 

SSU  Saybolt  Seconds  Universal 

t  time  (sec) 

T  temperature  (°F) 

T  reference  temperature  (°F) 

t  piston  thickness  (in) 

V  =  dV/dt  changes  in  volume  with  respect  to  time  (in  /sec) 

V  reference  volume  (in  ) 

o  2 

w  flow  area  coefficient  (in  /in) 

Xr  ram  position  (in) 

Xr  ram  velocity  (in/sec) 

2 

Xr  ram  acceleration  (in/sec  ) 

Xv  servovalve  spool  position  (in) 

Xy^  maximum  servovalve  spool  position  (in) 

a  cubical  expansion  coefficient  (°F~  ) 

8  bulk  modulus  (psi) 

3  liquid  bulk  modulus  (psi) 

X  absolute  viscosity  temperature   coefficient 

p  density  (slug/inJ) 

p  reference  density  (slug/in  ) 

pw50  density  of  water  at  60°F  (slug/in  ) 

^  absolute  viscosity  (lb  sec/in^) 

kinematic  viscosity  (in  /sec) 

wh  hydraulic  natural  frequency 
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I.  INTRODUCTION 

A.  GENERAL  DISCUSSION 

The  vast  majority  of  hydraulic  systems  currently  in  use 
are  of  a  dedicated  type.   That  is,  they  are  designed  for  a 
specific  function  with  consistent  and  well  defined  constraints 
In  most  cases  they  are  further  limited  to  a  loop  in  which 
the  pump  supplies  power  to  a  single  power  element.   Such 
systems  can  be  delivered  as  compact  units  with  a  minimum  of 
connecting  piping.   Some  examples  of  dedicated  systems  in 
Naval  use  are:   surface  ship  steering  units,  gun  mount 
train  and  elevation  units,  and  ammunition  conveyor  systems. 
Although  submarines  have  employed  dedicated  systems  in 
the  past,  all  modern  U.S.  submarines  use  constant  pressure 
systems  [Ref.  1,  Tab.  561-1].   The  constant  pressure  systems 
contrast  with  the  dedicated  type  in  that  the  supply  and 
return  headers  extend  the  entire  length  of  the  ship  and  serve 
a  multitude  of  power  elements.   The  total  piping  length  of 
these  systems  may  exceed  the  length  of  the  ship  by  a  factor 
of  four  or  more . 

B.  THE  NEED  FOR  A  LINE  MODEL 

For  the  packaged  configuration  of  the  dedicated  system, 
with  its  relatively  short  piping  lengths ,  the  current  practice 
of  neglecting  line  dynamics  and  lumping  the  system  volume 
has  produced  acceptable  analytical  results.   However,  even 
in  systems  with  low  line-volume  to  total-volume  ratios,  the 
computed  values  of  the  hydraulic  natural  frequency   (w^) 
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are  frequently  in  excess  of  40%  higher  than  measured  values 

LRef.  2  J.   Since  u>,  sets  the  upper  limit  on  system  response 
and  is  affected  by  the  volume  as  well  as  line  dynamics ,  a 
model  which  includes  these  factors  would  undoubtedly  produce 
more  accurate  results . 

It  was  the  intent  of  this  project  to  demonstrate  that 
the  submarine  hydraulic  system,  with  its  extremely  long  lines, 
requires  that  the  line  dynamics  be  included  to  achieve  any 
degree  of  accuracy  in  analytical  solutions. 

C.   POTENTIAL  FOR  SYSTEM  MODEL 

Discussions  between  the  author  and  Mare  Island  Naval 
Shipyard  Hydraulic  Design  Division  personnel  (code  260)  , 
as  well  as  a  Naval  Sea  Systems  Command  representative 
[Refs.  2  and  3]  ,  disclosed  that  there  is  currently  no  computer 
program  to  model  the  dynamic  behavior  of  a  submarine  hydraulic 
system.   Such  a  program  would  be  useful  in  the  design  of  new 
systems,  both  for  sizing  of  lines  and  elements  and  for 
determining  the  response  of  components  served.   For  existing 
systems,  it  would  allow  a  determination  of  the  effects  of 
alterations  prior  to  their  implementation. 

Additionally,  a  computer  program  of  a  dynamic  system 
model  would  be  advantageous  in  analyzing  vibrations  trans- 
mitted through  the  system.   The  "noise"  generated  by 
vibrations  is  a  critical  element  of  the  "Submarine  Sound 
Quieting  Program,"  while  the  fatigue  due  to  vibrations  is 
important  to  the  life  of  mechanical  elements  in  the  system, 
especially  the  seals  [Refs.  5,6,  and  7]. 
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The  development  of  a  dynamic  line  model  which  could  be 
programmed  in  a  modular  format  and  coupled  with  models  for 
the  individual  elements  would  facilitate  the  creation  of 
such  an  overall  system  model. 
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II.   SYSTEM  MODELED 

A.  MAGNITUDE 

It  was  recognized  early  in  the  undertaking  that  it  would 
not  be  possible  to  develop  a  model  for  an  entire  submarine 
hydraulic  system  in  the  time  available.   For  this  reason, 
the  stern  diving  planes  servo-loop  was  chosen  for  the 
primary  development  (Fig.  1) .   This  sub-system  was  considered 
to  be  representative  of  one  of  the  more  complex  loops  in  the 
overall  system.   (  a  loop  being  defined  as  the  fluid  path 
through  a  control  valve  -  to  the  power  element  -  back 
through  the  control  valve) .   Successful  modeling  of  the 
stern  planes  servo- system  would  readily  permit  expansion 
of  the  system  to  include  the  fairwater  planes  and  rudder 
loops  in  parallel  operations.   With  these  three  sub-systems 
modeled,   the  remainder  of  the  hydraulic  system  could  be 
relatively  easily  developed. 

B.  REALISM 

An  operational  ship's  hydraulic  system  was  deemed 
essential  as  the  basis  for  the  model,  as  it  would  allow  the 
computed  parameters  to  be  compared  with  actual  performance 
data.   The  hydraulic  plant  of  USS  PERMIT  (SSN  594)  was 
chosen  as  the  system  to  be  modeled.   It  was  considered 
representative  of  a  modern  submarine  hydraulic  system  and 
the  technical  information  was  readily  available  at  Mare 
Island  Naval  Shipyard. 
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III.      HYDRAULIC   FLUID 

A.  TYPE 

The  hydraulic  oil  used  in  all  U.S.  submarine  internal 
high  pressure  (3000  psi)  systems  is  Navy  symbol  2190TEP. 
This  is  a  petroleum  based  fluid  conforming  to  military 
specification  MIL-L-17331  [Ref.  8].   It  exhibits 
approximately  the  same  characteristics  as  an  SAE-30-Detergent 
motor  oil. 

B .  PROPERTIES 
1.   Density 

The  density  (p)  of  the  oil  was  obtained  from  Ref.  9 
in  terms  of  specific  gravity  (SG) .   The  values  supplied 
were 

SG  =  0.896     at  74°F. 

SG  =  0.0901    at  60°F. 

The  density  of  water  at  60  degrees  Fahrenheit  is 

-5         3 

9.3475  x  10   (slug/in  )  and  is  assumed  to  be  the  reference 

for  specific  gravity  unless  otherwise  specified  [Ref.  10] . 
Specific  gravity  is  defined  as  the  ration  of  the  density  of 
the  oil  at  the  given  temperature  to  that  of  water  at  60  degrees 
Fahrenheit. 

SG  -  P0/pW6Q  (1) 

Solving  for  the  density  of  oil  at  temperature  (T) 
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yields 

p  =   8.4221  xlO'5   slug/ in3   at  60°F 
p  =   8.3754  xlO"5   slug/ in3   at  74°F 

2.  Cubical  Expansion  Coefficient 

The  cubical  expansion  coefficient  (a)  is  the  fractional 
change  in  volume  of  the  fluid  due  to  a  change  in  temperature  and 
is  used  in  the  determination  of  density,  viscosity  (y,v), 
and  bulk  modulus  (S)  values. 

Density  values  arrived  at  in  the  previous  section 
were  used  to  obtain  a  value  for  the  cubical  expansion  coefficient 

Using  the  form  given  in  Ref.  11, 

p(P,T)  =  pQ[l  +  (P-PQ)/B  -  a(T-TQ)]  (2) 

where  the  subscript  (  )  denotes  reference  values  of  one 
atmosphere  of  pressure  and  60  F  temperature. 

Since  the  values  for  the  density  of  the  oil  were 
obtained  at  one  atmosphere 

equation  (2)  reduces  to 

p  =  Po  [1  +  a  (T  -  TQ)]  (2a) 

Inserting  appropriate  values  into  (2a)  and  solving  for  a 
yields 

a  =  3.9638  x  10"4  °F~1 

3.  Viscosity 

a.   Kinematic  Viscosity 

Reference  12  specifies  that  2190TEP  oil  shall 
have  a  kinematic  viscosity  (v)  of  190SSU  (Saybolt  Seconds 
Universal)  at  a  temperature  of  130  F. 
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Using  the  standard  conversion  [Ref.  14J  for 
kinematic  viscosity,  where  the  time  (t)  is  greater  than 
100  seconds  (SSU) 

v  =  0.220  t  -  (135/t)  (3) 

yields 

v  -  41.09  cs    at  130°F 

v  -  92.03  cs    at  100°F 

Converting  from  centistokes  (cs)  to  English  units  was 

accomplished  with  the  conversion  factor  (Ref.  15)  of  1.549  xlO 

2 

(in  /sec)/cs.   Thus  the  values  obtained  for  kinematic  viscosity 

were 

v  =  0.06367  in2/sec  at  130°F 

v  =  Q.1427  in2/sec  at  100°F 

b.   Absolute  Viscosity 

Knowing  values  of  density  and  kinematic  viscosity 
permits  the  determination  of  absolute  viscosity  (y) . 

y  =  pv  (4) 

Using  the  form  presented  in  Ref.  16, 

y  =  yQe  -  X  (T  -  TQ)  (5) 

at  one  atmosphere  of  pressure,  where  X  is  the  temperature 
coeffcient  of  absolute  viscosity  with  units  of  (  F  ) . 

Combining  equations  (4)  and  (5)   and  then 
solving  for  lambda  gives 

X  =  In  (P0v0/pv)/(T-TQ)  (6) 

Arbitrarily  setting 

T  =  100°F  and  T  =  130°F 
o 

in  conjunction  with  values  of  density,  and  kinematic 
viscosity  at  these  temperatures  yields 
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TQ  =100  T  =  130°F 

pQ  =  8.2886  xlO"5     p  =  8.1884  xlO"5 

vQ  =  1.427  xlO"1      v  =  6.367  xlO"2 
Substituting  these  values  into  equation  (6)  produced  a 
value  of 

X   =  2.731  xlO"2  °F~1 

Reference  17  presents  the  pressure  conversion 
for  absolute  viscosity  as 

log10(y/uQ)  =  cP  (7) 

However,  as  the  other  functions  are  related  by  the  difference 
in  pressure  divided  by  the  bulk  modulus  (AP/S) ,  it  would 
seem  reasonable  to  assume  that  c  in  equation  (7)  has  a  form 
of  (c-,/3).   This  would  then  yield  an  equation  for  absolute 
viscosity  at  constant  temperature  of 

]i  =yQe  -(ClAP/6)  (8) 

Lacking  any  other  information,  c-,  was  set  equal  to  one, 
and  a  final  form  for  absolute  viscosity  was  formulated  as 


where 


y  =  yQe((P/B)  -  A(T-TQ))  (9) 


p  =  p.p 
o 

P  =  reference  pressure  (o  psig) 

P  =  system  pressure  (psig) 

3  =  system  effective  bulk  modulus  (psi) 

T  =  system  temperature  (  F) 

T  =  reference  temperature  (  F) 
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4.   Bulk  Modulus 

Although  the  literature  on  hydraulic  systems  indicated 
that  the  bulk  modulus  (6)  would  be  supplied  with  the  oil, 
visual  inspection  of  shipping  containers  and  bills  of 
lading  proved  futile.   Likewise,  the  military  specifications 
(MIL-L-17331)  failed  to  produce  any  worthwhile  information 
in  this  regard.   A  nominal  value  of  bulk  modulus  for  the 
oil  was  therefore  assumed  to  be  220,000  psi  ,  as  recommended 
in  Ref.  18. 

C .   LIMITATIONS 

In  a  further  effort  to  maintain  the  scope  of  the  project 
within  reasonable  constraints,  it  was  decided  to  treat  the 
temperature  as  constant.   Inferences  drawn  from  the  Naval 
Ships'  Technical  Manuals  previously  cited  and  the  author's 
own  experience,   indicated  that  a  temperature  of  130  F  was 
a  reasonable  choice.   With  this  value  of  temperature,  the 
equations  for  the  other  properties  were  simplified. 

From  equation  (2) 

p  =  8.188  xlO"5  +  8.4221  xl0"5(P/f3)  (10) 

and  equation  (9)  yields 

y  =  5.213  xlO"6   rP/6)  .  (11) 
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IV.   MODEL  DEVELOPMENT 

A.   SYSTEM  OPERATION 

Submarine  steering  and  diving  control  surfaces  (rudder 
and  diving  planes)  are  positioned  through  mechanical  linkages 
by  hydraulic  rams.   Each  surface  has  normal  as  well  as 
alternate  sources  of  hydraulic  and  electrical  power  to  ensure 
reliability.   For  the  purpose  of  model  development,  only  the 
normal  mode  of  the  stern  planes  was  considered. 

The  normal  mode  of  operation  for  the  SSN  594  class 
of  submarine,  employs  a  position  control  electro-hydraulic 
servo-loop   (Fig.  1).  In  this  type  of  system  (Fig. 2),  control 
station  operator  manually  positions  a  stick  to  achieve  a 
desired  angle  of  inclination  of  the  planes.   This  produces 
an  electrical  signal  which  is  summed  with  the  feedback 
signal  from  the  position  output  to  create  an  error  signal. 
The  error  signal  is  amplified  and  transmitted  to  the  torque 
of  the  servovalve.   The  torque  motor  displaces  the  wand  spring 
of  the  servovalve-pilot ,  which  causes  an  imbalance  in  pilot 
pressure,  which  in  turn  shifts  the  main  spool  allowing  the 
servovalve  to  pass  oil  to  position  the  ram.   As  the  desired 
position  is  attained,  the  error  is  reduced  and  the  pilot 
returns  the  spool  to  a  neutral  positions ,  blocking  the  flow 
of  oil. 
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B.   SERVO- SYSTEM 

The  area  of  interest  for  the  model  development  was  the 
hydraulic  fluid  flow  in  the  servo-loop  as  depicted  in  Figure  3. 
By  inputting  the  valve  spool  motion  and  ram  load  force  as 
explicitly  defined  functions, the  influences  of  the  external 
factors  involved  in  the  feedback  loops  (the  servovalve  torque- 
motor  and  pilot  valve,  and  the  ram  loading)  were  minimized. 
This  permitted  the  generation  of  flow  and  pressure  equations 
around  the  loop  with  a  minimum  of  complications  due  to  the 
ancillary  characteristics. 

1 .   Servovalve 

The  two  servovalves  commonly  in  use  in  the  Submarine 
Force  are  both  manufactured  by  the  Sargent  Corporation  and 
are  designated  Sanders  Model  SV  438-10P  and  SV  438-15P 
[Ref.  19].   Design  and  performance  of  both  valves  are 
similar.   The  SV  438-10P  which  is  installed  in  the  (SSN-594) , 
was  used  for  this  simulation.   Figure  4  is  a  representative 
drawing  of  the  actual  valve  employed  -  a  four- land  four-port 
element.   However,  comparison  with  the  valve  depicted  in 
Figure  3  -  a  three-land  four-port  model  -  will  show  that 
the  flow  is  essentially  the  same  in  both  valves.   It  was 
considered  easier  to  follow  the  flow  paths  of  the  valve  in 
Figure  3 . 

All  technical  data  for  the  servovalve  was  obtained 
from  Ref.  20. 

a.   Flow  Paths 

The  pilot  valve  flow  was  disregarded, as  it  has 
its  own  supply  and  return  lines. 
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24 


at 
> 
i-i 
cfl 
> 

o 
> 

u 

cu 


cu 

•H 
fa 


25 


Designating  the  supply  flow  as  Q.    the  return 
flow  as  Q  ,  and  the  other  flows  by  their  respective  nodal 
points  (Fig. 3),  the  flow  equations  were  derived.   Using 
the  form  given  in  Ref.  21, 

E(^in  "  ZCW  =  dVo/dt  +  (V0/8)dP/dt         (12) 
the  results  were: 

Qin  =  Qx  +  Q2  (13) 

Qr  =  Q3  +  Q4  (14) 

Ql  =  Q3  +  Q5  +  Wsi  +  *i  <15> 

Q2  "  Q4  +  Q6  +  V2Ve2  +  ^2  •  <16> 

Here  the  volume  across  the  lands  of  the  valve  was  approximated 

as  zero.   The  superscript  (*)  represents  differentiation 

with  respect  to  time. 

Since  the  valve  is  manufactured  from  high 
strength  alloys  , 

vx  =  v2  =  0 

-     Vl  -  V2  • 
as  shown  in  Appendix  A. 

In  addition  to  summing  the  flows  through  the 
valve,  the  flow  rate  was  also  defined  through  the  lands  for 
three  different  valve  positions. 

(1)  Null  Positions.  In  the  centered  (null) 
position,  a  critical-center  servovalve  will  have  leakage 
past  the  lands.1   This  flow  is  described  in  terms  of  null 


!For  a  comprehensive  coverage  of  all  terms  used  in  this 
paper,  the  reader  is  referred  to  Ref.  2. 
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flow  pressure  coefficient  (K  ) ,  as 


CO' 

Q  =  %KcQ(AP).  (17) 

The  four  node  flows  for  X  =  0  were 

Ql  "  *Kco<Pin-p>  <18> 

^2  "  %Kco<Pin-p2>  <19> 

Q3  =  %Kco<PrPr>  (20) 

%   "  %Kco(P2-pr)  <21> 

with 

P3  =  P4  "  Pr- 
See  appendix  B  for  values  K   used  for  the  model. 

(2)   Valve  Spool  Off-Centered  Left.  The  valve 

travel  (X  )  was  arbitrarily  defined  as  positive  to  the  left 

(Fig. 3).   Employing  the  orifice  flow  equation  [  Ref .  22  ]  for 

the  fluid  passage  through  nodes  one  and  four,  and  an  empirical 

relationship  for  nodes  two  and  three,  the  flow  equations  were 

determined  to  be 

Qx  =  CdA/  2(Pin-P1)/p1  (22) 

Q2  =  %Kcoe-A(Xv/Xvm)2  (P.n-P2)  (23) 

Q3  =  V^W2  <Pl"Pr>  <24> 


Q4  =  CdA/2(p2-Pr)/Pr  (25) 

where  C,  =  0.61  for  full  periphery  rectangular  ports  and 

A  =  irdjX^  (26) 

with  d-.  defined  as  the  land  diameter  of  the  spool  valve. 

Since  a  value  of  leakage  flow  was  not 

determinable  from  Ref.  20,  the  empirical  relationships  of 

Appendices   C  and  D  were  developed  using  Ref.  23  as  a  guide. 
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(3)   Valve  Spool  Off -Centered  Right.  The 
equations  for  the  "negative"  flow  direction  were  similar  to 
those  developed  in  the  previous  section,  with  the  obvious 
exception  that  Q~  and  Q^  became  the  principal  flows.   In 
order  to  maintain  their  flows  as  positive  values,  X  was 
squared  and  included  in  the  radical. 

The  form  for  X  less  than  zero  was 

QX  -  %Kcoe-*(V3W2(Pin-pl>  <"> 


Q2  -  Cdw  /2X^(Pin-P2)/p2  (28) 


Q3  "  Cdw/2X^(P1-Pr)/p2  (29) 

Q4  -  %Kcoe-X(Xv/Xvm)2(P2-Pr)  (30) 


with 


W  =  7rd1  (31) 

and  all  other  variables  as  defined  in  previous  sections. 
2.   Lines 

The  line  equations  were  generated  with  the  intention 
of  accounting  for  all  factors  influencing  the  pressure  and 
flow  in  the  system.   The  schematic  diagram  of  the  actual 
system  employed  in  this  development  is  shown  in  Figure  5. 
For  development  purposes ,  only  the  odd  numbered  nodes  are 
shown  here  as  the  even  numbered  side  equations  were  identical, 
a.   Flow  And  Pressure 

Disregarding  all  but  the  main  line  piping,  the 
flow  path  from  the  servovalve  to  the  ram  was  as  depicted  in 
Figures  3  and  5  were  identified  by  the  same  numbers  .   The 
dynamics  of  these  segments  were  then  defined  at  each  point 
when  viewed  as  in  Fig.  6. 
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FIGURE  6  LINE  NODE  IDENTIFICATION 

Point  A  corresponds  to  nodal  point  (5)  of  Figure  3 
and  is  considered  to  be  located  at  the  entrance  of  the  pipe 
at  its  junction  with  the  servovalve.   Similarly,  point  B  is 
located  at  the  connection  of  the  two  different  diameter  pipes, 
in  the  entrance  of  the  downstream  pipe.   It  is  equivalent 
to  node  (7)  of  Figure  3.   Location  C  is  just  beyond  the  end 
of  the  line  segment  at  the  entrance  to  the  cylinder, 
identified  as  node  (9)  in  Figure  3. 

Utilizing  these  definitions  for  the  nodal  points 
and  viewing  the  direction  of  flow  as  always  being  from  A  to 
C,  the  dynamic  equations  for  each  element  were  developed. 

(1)   Node  Five .  The  pressure  and  flow  equations 
for  node  five  were  derived  from  Ref.  24. 

(32) 
(33) 


Q5  =  Qy  +  V57P/B57  +  V5? 
P5  =  Pl  "  P5K5°J/2A5 


The  double- subscripted  variables  were  employed 
to  denote  average  values  through  the  line  segment,  while  the 
single- subscripted  variables  represent  values  at  the  particular 
node .   Thus 
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V57  -  A5L5y  (34) 

with  L,-7  being  the  actual  length  of  pipe  from  node  (5) 
to  node  (7)  and 

A5  =  TTD5/4.  (35) 

As  can  be  seen  in  Figure  (6) ,  the  diameter  of  segment  A-B 
(5-7)  is  constant  and  equal  to  the  diameter  at  (5) .   The 
total  volume  and  actual  length  will  be  developed  further 
in  later  sections  of  this  report. 

The  average  pressure  in  the  interval  was 

P57  =  (P5  +  P7)/2  (36) 
with 

P57  =  /P57dt-  (37) 

The  change  in  volume  with  respect  to  time  (V) 

was  set  equal  to  zero  for  all  line  segments ,  because  the 

lines  are  made  from  high  strength  materials  which  do  not 

permit  any  significant  expansion  or  contraction.   The 

exception  to  this  would  be  in  cases  where  there  are  flexible 

sections  installed  (  Aeroquip  -type  hoses) .   There  are 

numerous  locations  in  the  system  where  these  are  installed, 

but  not  in  the  servo- loop  modeled.   In  the  instances  where 

flexible  sections  are  incorporated  in  the  system,  an  expression 

for  the  change  in  volume  would  have  to  be  developed  based  on 

the  properties  of  the  material  used. 

The  value  of  Kr  was  based  on  an  abrupt  change 

in  the  cross-sectional  area  of  a  pipe, and  the  development 

is  shown  in  Appendix  E.   As  used  in  the  above  equations,  the 

sign  of  Kc  changes  with  the  direction  of  flow  such  that  the 

second  term  of  equation  (33)  will  always  reflect  a  pressure 
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loss  for  the  referenced  point.   Thus,  when  flow  is  positive, 
as  shown  in  Figure  (6),  P,  will  be  greater  than  P,-,  but  when 
the  flow  is  reversed  (negative  in  Figure  (6))  Pr  will  be 
greater  than  P-,  . 

Q7  is  the  flow  at  node  (7)  and  will  be  defined 
in  a  subsequent  section.   P-,  is  the  pressure  in  the  main 
servovalve  cavity  and  is 

Px  =  JJjdt  (38) 

where  P-,   was  introduced  in  section  IV.B.l.a. 

(2)   Node  Seven.   The  flow  through  node  seven  is 
exactly  the  same  as  that  through  node  five,  with  subscripts 
changed  to  reflect  the  appropriate  line  segment. 

Q7  =  Q9  +  V79P79/B?9  +  V79  (39) 

P79  =  /p79dt' 

The  pressure  at  point  (7) ,  on  the  other  hand, 

must  reflect  the  losses  due  to  flow  in  line  segment  (57) . 

This  was  accomplished  with  a  combination  of  the  equations 

in  Ref.  24  with  those  in  Appendix  E:  (&m 

P?  =  P5  -  32u57Le57Q57/A5D5   -   1 .  Up^Q^/A^  -p7K?Q7/2A5  . 

The  second  term  on  the  right  side  of  equation 
(40)  is  the  Hagen-Poiseuille  term  for  fully  developed  laminar 
flow  in  pipes .   The  third  term  accounts  for  developing  laminar 
flow.   The  final  term  accounts  for  the  entry/exit  loss  at 
cross- sectional  area  changes,  as  shown  in  Appendix  E. 

All  three  terms  must  reflect  a  loss  of 
pressure  for  flow  in  either  direction.  Since  Q57  assumes 
a  negative  value  for  flow  reversal,  it  presents  no  problem. 
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2 
However,  Qcy  will  always  be  a  positive  quantity  unless  it 

is  defined  as 

q527  -Q5/|Q57I,  (41) 

in  which  case  the  sign  of  Qc7  will  reflect  the  direction 
of  flow.   As  noted  previously,  K-,   determines  the  sign  of 

the  final  term. 

2 
Ac  is  used  in  the  last  term  of  equation 

(40) ,  becaue  the  area  employed  is  always  that  of  the  smaller 

of  the  two  at  a  junction. 

Previously  undefined  variables  are 

Q57  =  (Q5  +  Q7)/2  (42) 

Lecy  =  effective  length  from  (5)  to  (7) . 

Effective  line  length  will  be  discussed  in  a  later  section. 
(3)   Node  Nine .   The  pressure  at  point  nine  has 

the  same  form  as  that  at  seven 

P9  =  P7-32y79Le79Q79/A7D^  -  1 . Up^Q^/A^-p^Q^A*     (43) 


with 


Q79  =  Q79" I Q79 I  <44> 


for  the  reasons  discussed  in  the  previous  section. 

Flow  through  this  junction  is 

Q-, ,  will  be  defined  with  respect  to  ram  flow  in  a  following 
segment  of  the  paper. 

Qel  =  CelP9  .  (46) 

The  coefficient  of  external  leakage  (C  -,)    is  the  allowable 
leak  rate  past  the  cylinder  tail  rods  to  the  atmosphere. 
Appendix  F  derives  the  values  used,  based  on  information 
contained  in  Ref.  25. 
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There  were  no  volume  terms  included  in 
equation  (45)  because  the  transition  from  the  end  of  the 
line  to  the  cylinder  entrance  was  considered  an  infinitesimal 
distance . 

b.   Secondary  Component  Effects 

It  is  readily  apparent  from  a  comparison  of 
Figure  (5)  with  Figure  (3)  that  the  simplified  system  omits 
numerous  valves,  fittings,  and  other  miscellaneous  elements 
which  are  installed  between  the  servovalve  and  ram.   Although 
these  "secondary  "  components  have  a  minor  influence  individually, 
when  combined  they  contribute  significantly  to  both  static 
and  dynamic  characteristics  of  the  system. 

The  approach  taken  in  this  development  was  to 
consider  them  in  three  categories:  branch  lines,  minor 
components,  and  major  elements. 

(1)  Branch  Lines .   Branch  lines  are  the  segments 
of  pipe  connected  to  the  main  line.   They  are  subjected  to 
system  pressure  but  do  not  ordinarily  allow  fluid  to  flow. 
Gage  lines,  hand  positioning  branches,  and  capped  sections 
were  included  in  Figure  (5) . 

These  lines  contribute  only  to  the  total 
system  volume. 

(2)  Minor  Components .    The  minor  components 
were  considered  to  be  any  elements  installed  in  the  system 
which  contributed  to  the  effective  length  and  volume  of 
the  system  but  were  not  specifically  listed  in  the  list  of 
major  components  [Ref.  26]. 
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The  elements  in  this  group  were  ball  valves, 
tees,  and  elbows.  Reference  27  was  used  as  the  source  for 
information  on  them  and  Appendix  G  explains  their  development 
for  use  in  the  overall  analysis  of  the  system. 

(3)   Major  Components.   The  deceleration  valves 
(decel)  and  the  blocking  valve  were  included  as  major  components 
These  were  listed  in  Ref .  26  with  specific  pressure  and  flow 
characteristics.   The  volumes  employed  were  a  "best  estimate" 
based  on  diameter  and  the  author's  personal  experience.   The 
effective  lengths  were  as  determined  from  the  procedure  in 
Appendix  H. 

The  function  of  the  blocking  valve  is  to 
transfer  to  an  alternate  source  of  hydraulic  supply  (the 
Vital  System),  in  the  event  of  a  loss  of  main  header  pressure. 
As  pressures  in  the  two  systems  are  ordinarily  the  same 
(within  lOOpsi.  of  each  other)  and  the  valve  has  internal 
seals,  cross-port  leakage  is  negligible. 

The  deceleration  valve  is  a  cam  operated 
full- flow  valve  which  is  installed  to  prevent  over- travel 
of  the  ram.   As  the  ram  approaches  its  maximum  limit  of 
travel,  a  collar  on  the  tail  rod  mechanically  trips  the 
decel  can,  shutting  the  deceleration  valve.   This  permits  a 
fluid  cushion  to  arrest  the  ram  travel  rather  than  allowing 
it  to  a  "hard- stop"  (metal  -  to  -  metal)  at  maximum  velocity. 
Motion  in  the  reverse  direction  is  permitted  by  flow  through 
the  internal  check  valve  until  the  linkage  allows  the  main 
deceleration  valve  to  reopen. 
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(4)   Combined  Effects.   The  various  elements 
were  combined  in  tabular  form  for  inclusion  in  the  system 
analysis  (Tables  1  and  2) .   These  were  input  directly  into 
the  computer  program  to  solve  for  the  various  parameters 
as  discussed  above. 

3.   Cylinder  and  Piston 

The  final  element  in  the  system  is  the  cylinder  and 
piston  combination.   Here  the  power  is  transmitted  from  the 
fluid  to  the  stern  planes  by  the  piston  and  its  associated 
linkages . 

a.   Fluid  Dynamics 

The  pressure  and  flow  equations  assumed  a  form 
similar  to  those  of  the  lines. 

(1)   Flow.   The  flow  in  the  cylinder  was  described 


by 


where 


with 


Qii  =  Q±£  +  *n  +  viipii/eii  <47> 

<kt   =  CU(P11-P12>  <48> 


Pll  =  /plldt  (49) 

Q. n  is  the  internal  leakage  from  the  high 

pressure  side  to  the  low  pressure  side  of  the  piston.   The 

determination  of  C. -  is  shown  in  Appendix  F.   Flow  was 

defined  as  positive  from  node  (11)  to  node  (12);  thus, when 

P-i  2  ^s  greater  than  P,,,  Q.   will  assume  a  negative  value 

as  required  by  the  line  equations . 
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vll  =  Allxr  (50) 

A-,-,  is  the  effective  piston  area 

All  =  ^ll  "  Dr)/4  (51) 

where  D,,  is  the  land  diameter  of  the  cylinder  and  D  is 

the  diameter  of  the  piston  tail-rod. 

r  is  the  velocity  of  the  ram. 


where 


Vn  =  AnLu  (52) 


Hi  "  Lo  +  Xr  <"> 


and  LQ  =  %(Lcy  -tp)  (54) 

L   is  the  length  of  the  cylinder  t 
is  the  thickness  of  the  piston.   Ram  position  was  defined 
as  positive  to  the  left  in  Figure  (5) ,  so  that  motion  to  the 
left  increases  V-,-,  while  movement  6f  the  piston  to  the  right 
decreases  V-,,  .   For  this  reason,  it  was  necessary  to  define 

V12  =  Lo  -  Xr  (55) 

and  ultimately 

*12  =  "Q«  -  *12  -   V12H2/612-  <56> 

Thus,  the  various  flows  in  the  system  would  automatically 

assume  the  proper  sign  to  reflect  the  direction  of  flow. 

(2)   Pressure.  The  pressure  equation  for  the 

cylinder  and  piston  arrangement  became 

PU  =  P9  -  32u11L11Qu/Aud21  -  1.14  P11«^1/a|1.  (57) 

There  was  no  entry  loss  term  because  it 

was  taken  into  account  with  Pq .   The  other  difference,  with 

regard  to  the  previous  line  loss  equations,  is  that  Qjiis  the 
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land  diameter  of  the  cylinder  while  A,-,  is  the  effective 
area  of  the  piston.   The  effective  area  was  used  because 
Q/A  is  a  velocity  term  employed  in  defining  viscous  loss. 
On  the  other  hand,  D-,,  was  used  to  account  for  the  surface 
area  upon  which  the  fluid  is  '"dragging";  i.e.,  the  walls  of 
the  cylinder.   Since  the  tail  rod  is  traveling  at  essentially 
the  same  velocity  as  the  fluid,  it  contributes  no  losses 
due  to  viscosity. 

As  in  previous  sections,  it  was  necessary 
to  define 

Qi!  =  Qii'lQnl  (58) 

to  account  for  flow  direction, 
b.   Piston  Dynamics 

Ram  motion  was  defined  in  terms  of  position  (X  ) , 
velocity  (X  ) ,  and  acceleration  (X  ) . 

Xr  =  /Xrdt  (59) 

Xr  =  JXrdt  (60) 

K   -  (AllPirA12P12AilVFL)/Mr     <61> 
The  load  force  (F,) ,  as  described  at  the  outset, 

was  input  as  a  ramp  function  dependent  on  ram  position  and 

velocity. 

The  flow  of  oil  (ni  .  -> )  term  was  included  to 

oil 

account  for  the  force  required  to  maintain  the  momentum  of 
the  oil.   It  was  determined  to  be 

"oil  =  p12*12      for  V°  (62) 

and 

m    =o   V         for  V°-  •        (63) 
moil   P11V11  r 
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The  total  mass  of  the  system  (Up)  was  defined 
as  the  mass  of  the  ram  plus  the  mass  of  the  oil  from  the 

low  pressure  side  of  the  piston  to  the  return  tank,  such  that 

tank 
Mr  =  mr  +  I  p±V.  (64) 

ram 

where  the  subscript  (^)  indicates  each  main  line  segment 
(flow  path)  from  the  ram  face  to  the  return  tank. 

C .   PROGRAM 

The  foregoing  equations  were  re-arranged  to  yield  the 
form  shown  in  the  attached  computer  program.   They  were 
modeled  in  a  Continuous  System  Modeling  Program  (CSMP  III) 
in  the  belief  that  it  was  the  best  modeling  method  available. 
However,  numerous  attempts  failed  to  accomplish  the  desired 
result. 

The  algebraic  equations  form  a  loop  for  which  no  solution 
was  found. 
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V.      SUMMARY 

A.  CONCLUSIONS 

It  is  felt  that  the  development  of  the  line  equations 
demonstrates  the  necessity  of  including  these  factors  in 
any  in-depth  analysis  of  the  submarine  hydraulic  system. 
It  is  also  the  author's  belief  that  the  approach  taken  in  this 
presentation  will  yield  worthwhile  results  for  any  system 
employing  extensive  quantities  of  piping. 

As  stated  in  the  introduction,  a  generalized  computer 
program  of  this  type  would  have  many  benefical  aspects  for 
both  the  submarine  designer  as  well  as  the  operator. 

B .  RECOMMENDATIONS 

The  first  recommendation  is  a  possible  means  of  clearing 
the  algebraic  loop  problem  currently  preventing  the  program 
from  executing.   This  would  require  clearing  all  of  the 
pressure  and  flow  equations  of  reference  to  any  variables 
except  for  themselves  or  ones  which  directly  invoke  an 
integral.   That  is 

P?  =  f(P79,P57,P1,P7,etc.)  . 
However,  it  could  not  take  the  form  where  P7  was  set  equal 
to  Pq  or  any  of  the  flow  terms  (Ch),  as  they  do  not  invoke 
an  integral  function. 

The  other  possibility  would  be  to  locate  another  program 
which  could  efficiently  handle  both  integral  as  well  as 
simultaneous  algebraic  equations. 
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APPENDIX  A 
Determination  of  Servovalve  Volume 

No  values  were  obtainable  for  the  volume  of  the 
servovalve.   Therefore,  the  following  method  was  employed 
to  obtain  a  realistic  approximation. 

Using  nominal  dimensions  given  in  the  technical  manual 
[Ref.  20  ] and  the  knowledge  that  the  valve  is  made  of 
aluminum  alloy,  a  volume  to  weight  ratio  was  determined. 

Dimensions    18.5  x  8.0625  x  9.5  in. 

Specific  Gravity  of  alloy  approximately  2.8 

Given  weight  69  lb. 

Density  of  water  at  60°F.   62.41b/ft 

Vbody  "  1417  in3 
Yalloy  =  0.101  lb/in3 

W  =  Vy  =  143.3  lb   for  a  solid  block 


then 


Wgiven/W  "  °-48 


Thus  the  valve  body  weighs   approximately  half  of  what  a 

solid  block  would  weigh.    It  was  therefore  assumed  that 

3 

the  volume  of  metal  in  the  valve  was  708in  .   Further, 

estimating  that  the  passages  for  the  pilot  system  occupy 

half  of  the  void  area  yielded  a  main  valve  volume  of  about 

3 

350  in  ,  and  since  this  volume  is  that  occupied  by  both  the 

supply  and  the  return  segments  of  the  valve,  half  was  assigned 
to  each. 

V1  =  V2  =  175  in3 


43 


APPENDIX  B 

Null  Flow  Pressure  Coefficient  Determination 
The  allowable  null  flow  leakage  was  given  as  600  cc/min. 

at  1200  psi.  [Ref.  20].   This  was  converted  to  standard  form 

3 
of  (in  /sec) /psi.  as  follows 

K   =  600 (cc.min)/ 1200  psi  x  1  in3/16.39cc  x  1  min/60  sec 

This  yielded 


and 


KcQ  =  5.0844  xlO"4  (in3/sec)/psi 


%KcQ  =  2.5422  xlO-4  (in3/sec) /psi 


APPENDIX  C 


Empirical  Valve  Leakage  Coefficient 
The  leakage  flow  across  a  valve  port  decreases  as 
the  valve  moves  to  block  the  port.   Since  a  plot  of  leakage 
flow  (Q-i)  versus  valve  position  (X  )  is  symmetrical,  with  an 
exponential  decay  as  shown  in  Fig.  C-l,  and  exact  data  was 
not  given  for  leakage  other  than  the  null  position  an 
empirical  expression  was  developed. 

'  Q£  MAX 


Figure  C-l  Leakage  Flow  vs  Valve  Position 
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The  most  appropriate  expression  to   fit   this   curve   took 
the   form  of 


2 
VQ£max=  e~UVXvmax) 


where 


Q„  =  leakage  flow  at  any  valve  position 

Q<?™„„=  maximum  leak  rate  =  %K 

x£max  a  co 

X„   =  maximum  valve  travel  from  null. 
Representative  values  were  then  chosen  such  that 

V<W  =  -02  at  VXvmax  "  -25 
That  is,  at  257o  of  maximum  valve,  travel  leakage  through 

the  "blocked"  port  of  the  valve  decreases  to  2%  of  the  null 

position  leakage. 

Solving  for  A 

X   =  ^ln)^max/QP>/<VXvmax>2 
X  =  62.6. 

In  reality  X  could  vary  from  near  zero  to  about  100. 

Realistically  it  would  have  a  value  between  50  and  100. 

X  =  65  was  chosen  as  the  value  to  employ  in  this 

presentation. 

The  form  used  was 


Q£ 


=  %K  e"A(Xv/Xvmax) 
co 


with  the  determination  of  maximum  valve  travel  (X  )  shown 

vm 

in  Appendix  D. 
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APPENDIX  D 


Determination  of  Maximum  Valve  Travel 


Reference  shows  that 


..-3 


at 


Qmax  =  80  GPM  =  80  x  231/60  iirVsec 


T  =  115°F 

APtotal  "  1200  Psi 
Ap =  600  psi 


port 
as  depicted  in  Figure  D-l. 


c 

,1 

S] 

L200  PSI 

I'l    1*1 

600  PSI 
Figure  D-l   Valve  Test  Flow  Path 

Using  the  orifice  equation 
Q  =  CdA(2Ap/p) 


% 


with 


Cd  =  0.61 
p  =  8.442  xlO 


-5 


A  =7Td-,X   where  d,  =  1  in. 

Xv  =  Xvm 


Q  =Qm 
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such  that 


Xvm  ■   (V*Cddl<2fiye>« 


resulting  in  X   =  0.0426  in. 

This  value  was  used  in  spite  of  the  fact  that  it  seemed 
large  for  a  ''standard0  pilot  controlled  servovalve. 
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APPENDIX  E 
Loss  Coefficient  Due  To  Abrupt  Area  Changes 

The  standard  form  of  the  pressure  loss  due  to  abrupt 
piping  diameter  change  is  given  as 

Hx  =  KV2/2g 
where 

H-,  is  head  loss,  expressed  in  feet  of  fluid 

Y  is  fluid  velocity 

g  is  acceleration  due  to  gravity  • 
Head  loss  can  also  be  expressed  as 
H1  =  Ap/a 

Y  =   Pg 
and 

Y  is  weight  density 
while  velocity  is 

V  =  Q/A  =  4Q/ttD2 
combining  the  above  equations 

AP  =  8PKQ2/tt2D4  . 
K  is  defined  from  Figure  D-l 
K  -  %(1-D1/D2)   for  +Q 
K  =  -(1- (D1/D2) 2)2  for  -Q. 
NOTE  D-.  is  always  the  smaller  of  the  two  diameters  being 
considered.  >-  Q+ 


D2         Dl 


Figure  E-l  Flow  at  an  Abrupt  Area  Change 
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APPENDIX  F 
Seal  Leakage  Coefficients 

Seal  leakage  was  based  on  permiss-able  values  listed  in 
Naval  Ship's  Technical  Manual, [Ref.  25.]  The  values  used  in 
this  presentation  were  for  •'used'"'  seals  with  maximum  allowable 
leakage  before  replacement  is  required.   This  was  felt  to 
be  reasonable  since  the  operating  temperature  of  the  oil 
is  normally  greater  than  the  test  temperature,  which  will 
produce  slightly  higher  lead  rates. 

The  leak  rates  were  given  as  one  milliliter  per  ten 
cycles  per  inch  of  rod  diameter  for  external  leakage,  and 
five  milliliter  per  five  minutes  per  inch  of  rod  diameter 
for  internal  leakage.   These  were  based  on  normal  operating 
pressures  and  a  fluid  temperature  of  100  F. 

Operating  pressures  were  assumed  to  be 

P  ,  =  2000  psi 

el        r 

?±l   =  1000  psi 

where  P  ^    is  the  pressure  for  external  leakage  and  P.n  is 

el        r  °  Li 

the  pressure  for  internal  leakage.   A  cycle  was  defined 
as  the  time  to  go  from  neutral  (zero  plane  angle)  to  full 
rise  to  full  dive  and  back  to  neutral  position.   For  the 
system  modeled,  stern  planes  travel  is  plus  or  minus  27°. 
Thus, a  full  cycle  would  be  four  times  this  value,  or  108°. 
The  normal  rate  of  travel  is  about  five  degrees  per  second. 
Thus,  108/5  =21.6  sec/cycle  or  approximately 
3  cycles/min. 
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Since  the  tail  rod  used,  had  a  diameter  of  four  inches, 
the  allowable  external  leakage  was  four  milliliters  per  ten 
cycles . 
Then 


lS£   x  3_cZ  =  l.2cc/min 


4ml  x 

lOcy  1ml   1  min 


and 


3 
Co  =  1.2  cc  x  1  min  x  1  in   x   1 

*       nan  ■-  60  sec  16.39cc    2000  psi 

-7     3 

Cel   =  6-1  xl°   (inJ/sec)/psi 

Similarly  the  internal  leakage  for  a  land  diameter  of 
8.25  inches,  rounded  to  nine  inches,  was  nine  cubic 
centimeters  per  minute. 
And 

C±l       60x16.39  x  1000 

CU  =  9'15   xl0~6    (in3/sec)/psi 
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APPENDIX  G 
Minor  Component  Analysis 

The  values  used  in  the  derivation  of  the  minor  component 
effects,  considered  herein,  were  obtained  from  Crane  Technical 
Paper  No.  410,  [  Ref.  27]. 

The  volume  was  determined  to  be 

V  =  AL 
where  A  was  the  cross-sectional  area  of  the  pipe  entering 
or  leaving  the  element  and  L  was  a  value  determined  for 
each  type  of  element.   The  values  for  L  were  chosen  as 

Ball  valves    L  =  3  x  D 

Tees  L  =  2  x  D 

Elbows         L  =  1.6  x  D 
In  cases  where  the  element  has  differing  entry  and  exit 
diameters,  D  is  the  larger  of  the  two. 

The  values  for  the  effective  lengths  were 

Ball  valves       L  -  3  x  D 

e 

Flow  through  tees  L  =  20  x  D 

°  e 

Branch  flow  tees  L  ■  60  x  D 

e 

45°  elbows        L  =  16  x  D 

e 

90°  elbows        L  =  30  x  D 

e 

where  the  diameters  were  as  above.   An  additional  factor  of 
K  (appendix  E)  would  have  to  be  included  for  diameter 
differences  in  the  main  flow  path. 

The  value  of  effective  length  cited  above,  were  determined 
for  Reynolds  numbers  (R  )  greater  than  1000.   Therefore, 
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it  was  necessary  to  modify  them 
Since         Re  =  4Q/(ttDv) 

and  L  ,=  L  R  /1000     for  R  <1000 

ea   ec  e  e 

with  ^ed=  t^ie  desired  effective  length 

L  =  the  calculated  effective  length  from  above 

then  Le  =  Lec  (4Q/ (ttDv)  )xl0~3 

If  R  > 1000  then   L  =  L 

e  -  e    ec 

so  that  for   4Q/ttDv>1000 

or  Qp/Dy>  785.4    L  =  L 

^K   -  e    ec 

and  for       Qp/Dy<785.4     L  =  1.273  xlO"3L  Qp/Dy 
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APPENDIX  H 


Effective  Length  of  Major  Valves 


The  standard  head  loss  equation 

o 

H-,  =  K  /2g  was  rearranged  in  the  form 

AP/pg  =  KQ2/2gA2 


This  yielded 


since 


K  =  2aPA2/pQ2 


A  =ttD2/4 

2    4     2 
K  =tt  APD78pCr 


and  since 


K  =  Le/D 
Le  =tt2APD5/8pQ2 
To  allow  Q  to  be  entered  in  the  form  presented  (GPM)  , 

a  conversion  was  included. 

3 
231  in  /gal  x  1  min/  60  sec 

so  that 

Q  (GPM)  x  231/60  =  Q  (in3/sec) 

and  Le  =  .0832  PD5/pQ2 
where 

D  =  diameter  of  valve  (in) 

P  =  given  pressure  loss  (psi) 

Q  =  given  flow  (GPM) 

3 
P  =  density  of  fluid  at  given  temperature  (slug/in  ) 
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APPENDIX  I 


TO  EMPLOY: 

SUPPLY  LINES;  Dl  IS  ALWAYS  THE  POINT  CLOSEST  TO  THE 
PUMP  -  THE  UPSTREAM  POINT 

RETURN  LINES;  Dl  IS  ALWAYS  THE  POINT  FARTHEST  FROM 
THE  TANK  -  THE  UPSTREAM  POINT 

DUAL  DIRECTION  (REVERSING)  LINES;  Dl  IS  THE  POINT 

NEAREST  THE  CONTROL  ELEMENT  (VALVE),  VIEWED  AS  IF 
THE  LINE  ALWAYS  SUPPLIES  PRESSURE  TO  THE  POWER 
ELEMENT 
ENDMAC 

INITIAL 

CONSTANT  KCCCC=600. ,FT=1200. ,LAMB=65. 

VALVE  LEAK  RATE  OF  MAIN  SPOOL 

NOTE:  LAMB  IS  STRICTLY  EMPIRICAL  AND  BASED  ON  AN 

AVERAGE  CURVE  OF  LEAK-RATE  VS  VALVE  DISPLACEMENT. 

CONSTANT  XVM=0 . Oij.25,  DLV=1 . 0 ,  CD=0 .  61 ,  Vl=175. ,  V2=175. 

CONSTANT  od5=i.90,th5=o.  25,  l5=i+8.o,ntb5=o.,ntt5=o.,  ... 

NIij5=0 . , NL95=1 . , NV5=0 . 
CONSTANT  0D6=1 . 9 ,  TH6=0 .  25 ,  L6=l+8 . 0 ,  NTB6=0 . ,  NTT6=0 . , . . . 

NLl;6=0 . ,  NL96=1 . ,  NV6=0 . 
CONSTANT   CD7=2 . 0 ,  TH7=0 .  27 ,  L7=228 . 0 ,  NTB7=2 . ,  NTT7=1 . ,  NLi+7=0 . , . . . 

NL97=6.,NV7=1. 
CONSTANT  0D8=2.0,TH8=0.27,L8=216.0,NTB8=3. ,NTT8=0. ,NLl+8=0. , . . . 

NL98=5.,NV8=1. 

CONSTANT  OD71=0.75,TH71=0.109,L71=8.0,NTB71=l.,NTT71=C.,... 

NLi|-71=0.  ,NL971=0.  ,NV71=0. 
CONSTANT  0D72=0 . 375, TH72=0 . 058 , L72=120 . 0 , NTB72=0 . , NTT72=0 . , . . . 

NLlj.72=0.  ,NL972=1.  ,NV72=1. 
CONSTANT  OD73=1.5,TH73=0.203,L73=12.0NTB73=0.,NTT73=0.,... 

NLi+7  3=0 . ,  NL97  3=0 . ,  NV7  3=0 . 

CONSTANT  OD8l=0.75,TH8l=0.109,L8l=8.0,NTB8l=l.,NTT8l=0.,... 

NLi;8l=0.  ,NL98l=0. ,  NV8l=0. 
CONSTANT  0D82=0 . 3750 , TH82=0 . 058 , l82=120 . 0 , NTB82=0 . , NTT82=0 . , . . . 

NLU82. ,NL982=2. ,NV82=1. 
CONSTANT  OD83=1.5,TH83=0.203,L83=12.0,NTB83=0.,NTT83=0.,... 

NLlj.8  3=0 . ,  NL98  3=0 . ,  NV8  3=0 . 

CONSTANT  PDR7=60.0,QDR7=50.0,RDR7=8.68E-05,VDV7=60. ,DDV7=1.09l+ 

CONSTANT  PBR7=50  .• ,  QBR7=50 . ,  RBR7=8 .  68E-05 ,  VBV7=20 . ,  DBV7=1 . i+6 

CONSTANT  E=3.0E07,POIS=0.3 

CONSTANT  B0=200000. 

CONSTANT  PR=50 . , PIN=3000 . , FA=0 . 1 , 0DRE=2 . 0 , THRE=0 . 1^8 , VR-2750 . 

CONSTANT  0ILCC=9 . 0 , CELCC=1 . 2 , PCEL=2000 . , PCIL=1000 . 

CONSTANT  DLR=8.255,DRR=lj..O,ODR=11.5,LCY=39.875,TR=5.125,MR=12.1 

CONSTANT   PI=3. 11+1593 
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KCC=0.5*KCOCC/(PT*983.1+) 

THE  VALUE  OF  (0.5)  IS  USED  BECAUSE  THE  TECH.  MANUAL 
LISTS  A  LEAK- RATE  FOR  THE  WHOLE  VALVE;   THE 
VALUE  SOUGHT  IS  THAT  FOR  HALF  THE  VALVE  (ONE  SET 
OF  PORTS). 

KCOCC  IS  TECH.  MANUAL  LEAK-RATE  IN  (CC/MIN). 

PT    IS  TEST  PRESSURE  USED  TO  ESTABLISH  THE  LEAK- RATE  (PSI), 

KCO   IS  LEAK-RATE  IN  ( (IN**3/SEC)/PSI) . 

AG=CD*PI*DLV*SQffTU. ) 

AG     IS  THE  PRODUCT  OF  THE  CONSTANT  COEFFICIENTS  USED  IN 
IN  THE  SERVO- VALVE  EQUATION. 

CFA-FA/1.1* 

fa    is  the  fraction  of  air  entrapped  in  the  system, 

expressed  as  a  percentage  of  the  total  system  volume. 

cfa   is  the  coefficient  of  the  air  term  in  bulk  modulus 
calculations. 

k5p,k5m=entco(di,d5) 
k6p , k6m=ent  co ( d2 , d6 ) 
k7p,k7m=entco(d5,d7) 
k8p ,k8m=entc0 ( d6 , d8 ) 

K9P , K9M=ENTC0 ( D7 , Dll ) 

k10p,k10m=ent  c0(d8,d12) 

these  coefficients  account  for  exit  losses  only: 
entry  losses  are  considered  as  exit  losses  for 
the  next  element  upstream. 

vp5,a5,lep5,d5la5,bcs5=vale(pi,e,pois,od5,th5l5,ntt5,ntb5,.. 

NLi|5,NL95,NV5) 
VP6 , A6 , LEP6 , D6LA6 , BCS6=VALE ( PI , E , P0 IS , 0D6 , TH6 , L6 , NTT6 , NTB6 , . 

NLi+6 ,  NL96 ,  NV6 ) 
VP7,A7,LEP7,D7,LA7,BCS7=VALE(PI,E,POIS,OD7,TH7,L7,NTT7,NTB7, 

NLi]7,NL97,NV7) 
VP8,A8,LEP8,D8,LA8,B(S8=VALE(Pl,E,POIS,OD8,TH8,L8,NTT8,NTB8, 

NLi+8,NL98,NV8l) 

V&B1,A7B1,LE7B1,D7B1,LA7B1,BCS71=VALE(PI,E,P0IS,0D71,TH71,.. 

L71,NTT71,NT37i,NLi|71,NL971,NV71) 
V7B2,A7B2,LE7B2,D7B2,LA7B2,BCS72=VALS(PI,E,P0IS,OD72,TH72,.. 

L72,NTT72,NTB72,NLl;72,NL972,NV72) 
V7B3,A7B3,LE7B3,D7B3,LA7B3,BCS73=VALE9PI,E,POIS,OD73,TH73,.. 

L73,NTT73,NTB73,NLl4-73,NL973,NV73) 

v8bi,a8bi,le8bi,d8bi,la8bi,bcs8i=vale(pi,e,pois,od8i,th8i,.. 

L8l,NTT8l,NTB8l,NLU8l,NL98l,NV8l) 
V8B2,A8B2,LE8B2,D8B2,LA8B2,BCS82=VALE9PI,E,P0IS,0D82,TH82,.. 

L82,NTT82,NTB82,NL^82,NL982,NV82) 


V8B3,A8B3,LE8B3,D8B3,LA8B3,BCS83=VALE(PI.P0IS,0D83,TH83,... 

L83,NTT83,NTB83,NLlj.83,NL983,NV83) 

B0I=1.0.B0 

BCII=2.0-"-(1.0+P0IS)/E 

BI1=B0I+BC1I 

BI2=BI1 

B15=BI1 

BI6=BI1 

BI7=BI1 

Bl8=BIl 

BI9=BI1 

BI10=BI1 

BI11=BI1  ■ 

BI12=BI1 

BCDV8=BI1*VDV8 

BCDV7=BI1*VDV7 
BCBV8=BI1*VBV8 

BCBV7=BI1*VBV7 

Bl57=BOI+BCS5/V57 
BI68=B0I+BCS6/V68 

BI79=B0I+  ( BCS7+BCS71+BCS72+BCS73+BCS714-+BCDV7+BCBV7 )  /V79 
Bl8l=B0I+(BCS8+BCS8l+BCS82+BCS83+BCS8^+BCDV8+BGBV6)/V-::-l 

cq57,cq57s=linco(a5,d5,Od5) 
cq58 , cq68s=linc0 ( a6 , d6 , 0d6 ) 

CQ79 , CQ79S=LINCO ( A7 , D7 , 0D7 ) 
CQ81 , CQ8lS=LINC0 ( A8 , d8 , 0D8 ) 
CQ11,CQ11S=LING0(A11,D11,0DR) 
CQ12 , CQ12S=LINC0 ( A12 , D12 , ODR ) 

PlZ=(PIN-PR)/2,0 

P2Z=01Z 

P57Z=P1Z 

P68z=P1Z 

P79Z=P1Z 

P81Z=P1Z 

P11Z=P1Z 

P12Z=P1Z 

Q7Z=0.0 
Q8z=0.0 

Q9Z=0  .0 

Q10Z=0.0 

Q11Z=0.0 

Q12Z=0.0 

Q57Z=0.0 

Q68z=0.0 

Q79Z=0.0 

Q8lZ=0.0 

All=PI*(Dll*Dll-DllR*DllR)A.O 
A12=PI-*(Dl2-::-Dl2-D12R*Dl2R)A.O 
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D11=DLR 
D12=DLR 
D11R=DRR 
D12R=DRR 

a5s=a5-»-a5 

A6S=A6*A6 

A7S=A7*A7 

a8s=a8*a8 

L57=LA5 
A68=LA6 

L79=LA7 
L8l=LA8 

lep57=lep5 
lep79=lep7+ldv7+lbv7 

LEP68=LEP6 
LEP8l=LEP8+LDV8+LBV8 

v57=vp5 
V68=VP6 

V79=VP7+V7B1+V7B2+V7B3+VDV7+VBV7 
V81=VP8+V8B1+V8B2+V8B3+VDV8+VBV8 

LDV7=SPECV ( PDR7 , QDR7 , DDV7 ) 

LDV8=LDV7 

LBV7=SPEGV ( PBR7 , RER7 , DBV7 ) 

LBV8=LBV7 

VBV8=VBV7 

VDV8=VDV7 

DRE=0DRE-2.0-::-THRE 

BIR=B0I+2.0-::-((l.+P0IS)-«-0DRE-"-0DRE+(l.-P0IS)v-DRE-::-DRE)/... 
(E*(ODRE*ODRE-ORE*ORE) ) 

CEL=CELCC  /(PCEL*983.1j.) 
CIL=CILCC/(PCIL-"-983.i;) 

XRZ=0.0 
ZRDZ=0.0 

Dl=3.-«-DLV 
D2=D1 

LC=0.5-»-(LGY-TR) 

DYNAMIC 
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*#   PROCEDURE  SERVO-VALVE  **# 

PROCED  Q1,Q2,Q3,QI|-=SERV0(AG,KC0,LAMB,PIN,PR,P1,P2,... 

XV,XVM,RH0l,RH02,RHCR) 

CEX=KCO*EXP  ( -LAMB*  ( XV/XVM) **2 ) 

IF(XV)10,20,30 
10    Ql=CEX*( PIN-PI) 

Q2=AG*SQRT  ( XV*XV*  ( PIN-P2 )  /RH02 ) 

Q3=AG*SQRT ( ZV*ZV* ( P1=PR ) /RHOR ) 

Ql|=CEX*(P2-PR) ' 

GO  TC  kO 
20   Ql=CEX*(PIN-Pl) 

Q2=CEX*(PIN-P2) 

Q3=CEX*(Pl-PR) 

Ql|=CEX*(P2-PR) 

GO  TO  kO 
30   Ql=AG*XV-*SQRT  (( PIN-PI )  /RHOl ) 

Q2=CEX-"-(PIN-P2) 

Q3=CEX-"-(  PI-PR) 

Q[|=AG(XV(SQRT(  (P2-PR)/RH0R) 
1|0    CONTINUE 
ENDPRO 

#*   PROCEDURE  KPRO  «*» 

PROCEED  K5,K6,K7,K8,K9,K10=KPR0(P1,P2,P5,P6,P7,P8,P9,P10,K5P,... 
K5M,K6P,K6M,K7P,K7M,K8p,K8M,K9P,K9M,K10P,K10M) 
if(pi.ge.p5)go  TO  100 
k5=k5m 

go  TO  101 

ioo  k5=k5p 

101  IF(P2.GE.P6)G0  TO  102 
K6=K6M 

GO  TO  103 

102  K6=K6p 

103  IF  (P5.GE.P7)GO  to  1014. 
K7=K7M 

GO  TO  105 
101+  K7=K7P 

105  if(p6.ge.p8)go  to  106 
k=k8m 

GO  TO  107 

106  K8=K8P 

107  IF(P7.GE.P9)G0  TO  108 

K9=K9M 
GO  TO  109 

108  K9=K9P 

109  IF(P8.GE.P10)G0  TO  110 

K10=K10M 
GO  TO  111. 

110  K10=K10P 

111  CONTINUE 

ENDPRO 
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***   PROCEDURE  EFFECTIVE  LINE  LENGTH  «*» 


PROCED  LE57=PR057 ( LEP57 , L57 , Q57 , D5 , RH057 , MU57 ) 
CLE57=Q57*RHOS7/ (mu57*d5 ) 
IF(CLE57.LT.785.i;)GO  TO  200 
LE57=L57+LEP57  " 
GO  TO  201 

200  LE57=L57+1 . 27  3E-0  3*CLE57*LEP57 

201  CONTINUE 
ENDPRO 

PRO CED  LE68=PR068 ( LEP68 , L68 , Q68 , D6RH068 , MU68 ) 
CLE68=Q68RH068/ (MU68*d6 ) 
IF(CLE68.LT.785-1;)G0  TO  202 
LE68=L68+LEP68 
GO  TO  203 

202  LE68=L68+1.273E-03*CLE68*LEP68 

203  CONTINUE 
ENDPRO 

PROCED  LE79=PR079(LEP79,L79,Q79,D7,RH079,MU79) 
CLE7  9=Q7  9-"-RH07  9/  ( MU7  9-*D7 ) 
IF(CLE79.LT.785.U)GO  TO  20k 
•  LE79=L79+LEP79 
GO  TO  205 
20k  LE79=L79+1.273E-03*CLE79*LEP79 

205  CONTINUE 
ENDPRO 

PROCED  LE8l=PR08l(LEP8l,L8l,Q8l,D8,RH08l,MU8l) 
CLE8l=Q8l-::-RH08l/  (MU8l*D8 ) 
IF(CLE8l.LT785.1|.)G0  TO  206 
LE8l=L8l+LEP8l 
GO  TO  207 

206  LE8l=L8l+1.273E-03-»-CLE8l*LEP8l 

207  CONTINUE 
ENDPRO 

***  PROCEDURE  MASS  OF  OIL  *-"-"- 

PROCED  MD0IL,MT=M0PR0(MR,XRD,RH011,RH012,V11,V12,V11D, 

V12D,RH057,RH068,RH079,RE081,RH01,RH02,V1,V2,V57,V63, 

V79,V81,... 

RHOR,VR) 

IF(XRD.GE.O.)GO   TO   ij.00 
K0IL=RH0ll*Vll+RH079J.:-V79+RH057*v57+RH0*Vl+RH0R-::-VR 

GO  TO   i±01 
ij.00  M0IL=RH0l2#Vl2+RH08l*v8l+RH068#RH02*-V2+RH0R*VR 

MD0IL=RH0l2v-V12D 
1+01  MT=MR+KOIL 

ENDPRO 
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P1D=(Q1-Q3-Q5)/V1*B1I) 
P2D=(Q2-Q^-Q5)/(V2*B2I) 

P57D=(q£-Q7)/(b57I*v57) 
P68d=(Q6-q8)/(b68i*v68) 

P79D=(Q7-Q9)/(B79I*V79) 
P81D=(q8-Q10)/(b81I-::-v81) 
PllD=(Qll-OlL-VllD)/(VllI-rrVll) 
P12D=(Q12+QIL+V12D)/(B12I*V12) 

P5=2.0^P57-P7 

P7=2.0*P79-P9 

P9=(Q9-QH)/CEL 

P6=2.0*p68-p8 

P8=2.0*p8l-P10 

P10=(Q10-Q12)/CEL 

«*#  line  and  ram  flows  «h* 
q5=a5*sqrt(2.o*(pi-p5)/(RHo5*k5)) 

q6=a6*SQRT ( 2 . 0* ( P2-P6 ) /  ( RH06*K6 ) ) 

THE  FOLLOWING  EQUATIONS  IN  IMPLICIT  FORM  WERE  INSERTED  IN 
AN  EFFORT  TO  CLEAR  THE  ALGEBRAIC  LOOP  PROBLEM.   IT  DID  NOT 
WORK. 

THE  NEXT  LOGICAL  STEP  WOULD  BE  TO  COMPACT  ALL  EQUATIONS 
SUCH  THAT  EACH  VARIABLE  IS  EQUATED  ONLY  TO  ITSELF  OR  A 
VARIABLE  WHICH  IS  DIRECTLY  COUPLED  TO  AN  INTEGRAL. 

Q7=IMPL(Q7Z,0.001,Q7CC) 

Q7C=(  (p5-P7-CQ57S-"-RH057-"-ABS(Q57(-0.5-«-RH07v-K7Q7Q7/Q5s)  . . . 

/(CQ57*mu57*LE57))*2.-q5 
Q7CC=Q7C 

q8=impl(q8z,o.ooi,q8cc) 

Q8c=(  (P6-p8-CQ68s-"-RH068-::-q68-u-ABS(q68)-0.5-"-RH08wK8-::-q8/a6s)  . . . 

/  ( CQ68#MU68*-LE68 ) )  *2 .  -Q6 
Q8CC=Q8C 

Q9=IMPL(Q9Z,0.001,Q9CC) 

Q9C=(  (P7-P9-CQ79S-::-RH079-»-Q79-~-AES(Q79)-0.5-"-RH09-:'-K9-«-Q9-»-Q9/A7S. . . 

.  '  / ( CQ79*MU79*LE79 ) )*2 . -Q7 
Q9CC=Q9C 

Q10=IMPL(Q10Z,0.001,Q10CC) 

Q10C=  ( ( P8-P10- CQ3lS-::-RH08l-::-Q8l-:t-AES  ( Q8l ) -0 . 5*RHOlO-::-KlO*QlO*QlO . .  . 

/A8s)/(CQ8l-::-MU8l*LE8l)  )*2.-Q8 
Q10CC=Q10C 

Q11=IMPL(Q11Z,0.001,Q11C) 

Q11C=  ( P9-P11-CQ11S-::-RH011-::-Q11v-ABS ( Qll ) ) /  ( CQll*KUll*Lll ) 

Q12=IMPL(Q127,0.001,Q12C) 

Q12C=  ( P10-P12- CQ125«-RH012«-Q12*ABS  ( Q12 ) ) /  ( CQl2-:c-MUl2-::-Ll2 ) 
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qS7C    /(cqS7*mu57*leS7) 
q57cc=qS7c 

g|t^^!w68»Kli§68S^((l68)-0.5»BH08«8»Q8«i6/A68)... 

/(CQ68*KiU68*LE68) 

q68cc=q68c 

^C~(P7-pK  ' '  • 

Q7  /(CQ79*KU79*LE79) 

Q79CC=Q79C 

Q51C-tP0  ^s))(cq8i..MU8i-le8i) 
q8icc=q8ic 

QIL=CIL*(PH-P12) 

b1i=cfa/p1+bi1 
b2i=cfa/p2+bi2 
bri=cfa/pr+bir 
b5i=cfa/p5+bi5 

B6l=CFA/P6+BIb 
B7I=CFA/P7+BI7 
B8l=CFA/P8+Bl8 
B91=CFA/P9+BI9 
B10I=CFA/P10+3I10 

BllI=CFA/PH+BIll 
B12I=CFA/P12+BI12 
b57I=CFA/p57+BI57 
B68l=CFA/P68+BIb8 
B79I=CFA/P79+BI79 
B8lI=CFA/P8l+Bl8l 

****DENSITY  VALUES  «H» 
pwop=8..256E-05+8.^23E-05-PR-»-BRI 

RH06=8   256E-0$+8.i|23E-05-"-P6-::-B6l 

ptoM.256e-0$+8.U23E-05*p8»B6I 

RH09=8   256E-05+8 . fi23B-05»P8»B9I 

RH011=8 .256E-05+8.^23E-05*Pll*BllI 
RH012=8  256e-05+8.1l23E-0^-P12-"-B12I 
RH057=8   2^E-05+8.U23E-05-::-P57-::-B57I 

BH079-8  2565-01+8.1*2335-0^7  *B7?I 

RH081=8 .  256E-05+& .  l|23E-05»P8l»B8lI 


_ 


***  VISCOSITY  EQUATIONS  *** 

MU57=8 .  198E-06*EXP  (  P57*B57I ) 
MU68=8 . 198E-06*EXP ( p68*B68l ) 
MU79=8. 198E-06-::-EXP  ( P79-"-B79I ) 
MU81=8.198e-06*EXP(p81*b81I) 
MU11=8 .  198e-06*EXP  ( P11---B11I ) 
MU12=8 .  198E-06-::-EXP  ( P12*B12I ) 
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